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ABSTRACT 

The possible conversion of neutron matter to strange matter and its relevance to 
astrophysical problems are discussed. Particular attention is paid to the form in which 
this hypothetical "burning" process should be propagated in dense matter and to the 
moment when it should start. 
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1. INTRODUCTION 



Since first stated by Witten [1] the hypothesis of strange matter (hereafter SQM) being 
the true ground state of cold hadronic matter has attracted considerable attention [2]. 
Among the several aspects studied within the strange matter hypothesis we may mention 
the problem of the conversion of (cold) neutron matter to strange matter in astrophysical 
objects. This has been first studied in [3,4] and refined in [5,6] and all these works have 
assumed that the process is analogous to an ordinary (slow) combustion. This is of course 
very reasonable in a first approximation, but unfortunately things seem not to be that 
simple. Not only the actual way in which the conversion develops is subject to considerable 
doubt, but the moment when it can start is also controversial. This point turns out to 
be very important because it can alter dramatically the fate of the host star. We shall 
attempt to summarize in this work what is known (and unknown) about this problem, 
which appears to us far from being closed. 

2. HOW AND WHEN DOES THE BURNING START ? 

A variety of mechanisms for the conversion to start have been discussed in [7] , tentatively 
classified as "primary" (in which SQM is produced somehow inside the dense environment) 
and "secondary" (SQM or a large amount of energy to produce it are seeded from the 
outside of the star). Clearly, the secondary mechanisms, even though perfectly viable, are 
much difficult to assess and and we shall not discuss them here. 

Perhaps the most likely primary mechanism discussed in the literature is the 
conversion via two-fiavor quark matter. The scenario relies on the possibility of deconfining 
the light quarks from ordinary nucleons as the stellar matter is compressed/cooled in 
a supernova aftermath. This is sketched in Fig.l ; the transition point is achieved at 
a density ric and the quark region relaxes to the SQM state on a timescale typical of 
weak interaction r-uj lowering its energy per baryon number unit. From the "macroscopic" 
point of view, and since Ty, <^ Tdynamical the process may be modelled as isocore. The 
released energy difference builds up an overpressure AP whose whose value can be found 
for a given model calculation. In [8] we have employed a free massless quark gas + bag 
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pressure and the Bethe- Johnson I [9] equation of state for neutron matter to calculate 
AP ~ 70ll0 MeV fm~^, depending somewhat on the parameters. This result shows 
that the overpressure is huge and its decay may be very important dynamically (see below). 
Note that, if at all, this process must occurr shortly after the gravitational collapse itself 
since a significative compression is achieved on a timescale r ~ 1 s [10]. 



Fig.l. A possible route for SQM appearance viewed in the free-energy vs. density plane. 
Deconfinement is achieved at a value ric (dashed line) and relaxation to SQM is further 
achieved on a weak-interaction timescale (pictured with a wavy line) . 

If SQM is not produced as the result of two-flavor quark matter "relaxation", it is 
possible that it could arise because of thermal nucleation of SQM bubbles inside hot nuclear 
matter. The idea is that if a single bubble of the fluctuation distribution grows beyond a 
critical size [11] it will convert the sourroundings afterwards. A simple calculation using a 
kinetic Fokker-Planck-Zel'dovich approach has been presented in [12]. The growth of the 
bubbles has been described by the expression 



dt dr ^ ' 



where /(r, t) is the time-dependent size distribution of bubbles and ^ is the nucleation 
rate. The appropiate solution for this problem may be written as 
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e = 2.2 10-2 {T/a)^ n}c t-^ exp{-3.1 N^c a/T) 



(2) 



with Nqc the number of quarks of the critical size bubble (< 300), a weak- 
interaction timescale and a the surface tension coeflBcient. A lower bound to the 
temperature at which at least one bubble larger than the critical size will nucleate can 
be found by imposing 



where At is the time available for prompt nucleation (~ 1 s) and Vnuc ~ (1 km)^ 
is the volume in which the nucleation can take place. As a < (70 MeV)^ from detailed 
calculations [13] and fittings to phenomenological models [14], the condition (3) can be 
re-written equivalently as 



This bound can be certainly satisfied in the first seconds of the neutron star life but 
not after that since the initial cooling phase is dominated by neutrino emission, which is 
in turn a strong function of T. Therefore, we again conclude that the conversion is likely 
to start immediately after the gravitational collapse itself. The same conclusion holds if 
the conversion is triggered by "dormant" strangelets [15] which become active as soon as 
a full neutronization of the core is possible (r ~ Is). 

3. WHICH IS THE PREFERRED COMBUSTION MODE ? 

As we have discussed in section 2, a short timescale ~ 1 s for SQM to be "operative" 
inside a proto-neutron star appears to be a general feature of any triggering mechanism. 
This is, of course, not enough to address the aftermath, since the form in which the 
conversion front propagates remains unclear. For a given chemical potential diff"erence A/x, 
the conversion n — > SQM releases energy much in the same way as a chemical reaction 



^ At Vnuc > 1 



(3) 



T > 2 MeV 



(4) 
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does. Therefore, the well-known theory of combustions [16] can be applied to study the 
process. Physically, steady state combustions can occur in two modes : deflagrations and 
detonations. Deflagrations (or slow combustions) are subsonic burnings in which quasi- 
equilibrium among the burnt and unburnt fluids holds. Detonations (or rapzd combustions) 
are supersonic and can be described as a shock wave followed by the combustion front. 

Assuming that a slow combustion occurrs, the velocity of the flame has been calculated 
to be Vs ~ few km/s at most [3,4]. These works have focused mainly on the reaction 
kinetics of the conversion, but the actual boundary conditions in a star and the gas 
motion should play a role as well. To study these effects, a linear stability analysis of 
this propagation mode has been performed [17] in which the stability condition was found 
to be 



with j the mass flux onto the flame, a the surface tension, g the local gravitational 
acceleration and pi,p2 the fluid densities of the unburnt and burnt sides respectively. 
On the other hand, by deflnition [16] 



from which either P2 > Pi and p2 > pi (detonations) or P2 < Pi and p2 < pi 
(deflagrations). The important point is that if we impose the latter conditions [18], then 
eq.(5) can never be satisfled. We conclude that deflagrations are always Rayleigh-Taylor 
unstable in dense matter for any perturbation wavenumber at the linear level. 

These results indicate that, unless the instability can be controlled in some way (e.g. 
by non-linear effects), the flame should become turbulent and accelerate. Laboratory 
observations show that, in such a case, the foldings of the front increase, thus increasing 
the area and the local burning rate. The heat transfer mechanism switches from thermal 
conduction to turbulent convection, which gives a much larger velocity Vg- An estimate of 
the latter can be given in terms of a lenght scale, taken to be the radius of the spherical 



/ < 4:agplpl 
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(5) 



(P2 - Pi) 




{P2 - Pl)piP2 
{P2 - Pi) 



(6) 
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front r (the only available lenght of the problem) divided by the timescale for the instability 
to grow trt- Given that trt niust necessarily be some multiple of the dynamical time 
(G'p)-i/2 we find 



r 



10^ ( 



r 



) cm s 



-1 



(7) 



'^turb — 



1 km 



which probably contains an uncertainty of a factor of 10 in either direction. The point here 
is that not only Vturb ^ ''-'s but also it is likely to be comparable to the speed of sound Cg at 
some finite radius. However, since the turbulent convection burning must also be subsonic, 
its description is also contained in the lower branch of the Hugoniot adiabatic. Note also 
that, due to the supranuclear density values of the stellar matter, trt — 1 ms and thus 
the acceleration process will be very fast. Moreover, a further change of the combustion 
mode is to be expected, as we shall now discuss. 

We have previously noted that boundary conditions must be taken into account for a 
reliable and consistent description of the fiow. As is well-known [16], these boundary 
conditions (analogous to the "closed pipe" ones) force the existence of a shock wave 
propagating ahead of the flame. The strenght of the shock is uniquely determined by 
the velocity of the flame and the state of the unburnt material. Now, if the flame 
accelerates the shock becomes stronger, which in turn means that the neutrons are more 
and more pre-heated and pre-compressed before meeting the flame. Such pre-heating 
and pre-compression effects on the matter entering the combustion zone help to further 
accelerate the burning rate and establishes a positive feedback loop [19]. 

The crucial point is therefore to establish which is the velocity of the flame for which 
the shock is strong enough to "ignite" the neutron matter by itself. This type of calculations 
have been longly performed for laboratory gases (for planar and spherical geometries) and 
it has been repeatdely found [19] that for vturb — 0.1 Cg the so-called auto-ignition 
condition is met and thus the only possible propagation is the one corresponding to the 
Chapman-Jouget point on the Hugoniot adiabat [16]. 

Very recent work [20] on this problem has revealed another important feature of the 
n — > SQM combustion : it has been found that deflagrations are not allowed unless 
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Fig.2. a) Laminar initial stage of the deflagration, arrows indicate the velocity of the flame 
and pre-compression shock b) folding of the flame due to instability grow c) fully developed 
turbulent flame (note the acceleration of the pre-compression shock —broken arrow—) 

Pi < 1.63 po ( Po being the nuclear saturation density). The reason for this is simply 
that for higher densities the velocity of the flame becomes larger than the velocity of the 
pre-compression shock. Physically this means that for this regime the flame and the shock 
coalesce in a single region, a feature that characterizes a detonation wave [19]. However, 
given the short timescale for the flame to become turbulent, it is very likely that the 
detonation develops a width S ~ r^Cs < 1 m. So, even though the kind of combustion 
is limited by the the weak interactions rate (producing s quarks), the heat transport 
mechanism must be recognized as a fundamental ingredient in the problem. It appears 
from the above results that we should consider ab initio a (turbulent) detonation mode. 

To conclude this section two remarks are in order. First, some doubt has been 
expressed about the ability of "slow" weak interactions to sustain an energy release fast 
enough for a detonation to propagate. We note that, even for the longest reasonable 
timescale of ~ 10~® s the front thickness is always negligible with respect to the dimension 
of the burning region (~ km) and, much more important, with respect to the lenghtscale 
of significative variations of the pressure, etc. estimated as L = [dlnP/dr)) ^ . In fact 
8 -C L is a sufficient condition for the stability of a propagating detonation front [21], 
or in other words, "slow" weak reactions are fast enough to sustain a detonation unless 
(say) L < 10 — 100 m . To be more quantitative, however, a complete numerical 
simulation, including the (strong) dependence of the weak interactions rate with T [6] 
should be carefully performed. The second remark is that it is by no means necessary 
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that the burning passes through an initial slow combustion stage and then accelerates. 
On the contrary, a "direct" initiation is more likely if SQM is produced by the decay of 
a two-flavor plasma in bulk (section 2). In this case, the overpressure AP decays into a 
shock wave with strenght 5P c^l AP/2, which should be enough to ignite the matter by 
itself (see [22] for a through discussion of this and related points). 

4. CONCLUSIONS 

We have discussed in this talk several topics related to the possible appearance and 
propagation of SQM inside dense neutron matter. We have argued that 

a) SQM is expected to form/begin to "operate" after ~ 1 s following the collapse at most, 
and we are not aware of any mechanism to postpone the conversion until, say, ages ~ yr 
(so that all neutron stars must be strange stars [1,2,7,12,18]). 

b) Even if a deflagration may develope first [3,4], it should become turbulent and switch 
to a detonation shortly after. Alternatively, a detonation may be directly initiated. 

The conclusions (to which a factor < / < 1 gauging the uncertainties of the 
problem may be added) point to the need for a better understanding of compact star 
interiors, specially for the role of few-quark states. Preliminary steps have been taken 
in [23]. It is also very important to improve the schematic model in which the sudden 
energy release from a detonation affects the supernovae events (hopefully in a positive way 
!) [8,24]. The conceptual similarity of the latter with the "delayed detonation model" for 
type I supernovae [21] is quite remarkable and may merit a closer look. Note, however, that 
even if a detonation is not formed after all, the very presence of a strong turbulent flame 
propagating at high (but subsonic) velocity may be more than enough to add substantial 
energy to a stalled prompt shock. 
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